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Abstract. It is now about 10 years since the evidence, based on Type la supernovae, for the 
acceleration of the expansion of the Universe was discovered. I will discuss some aspects of the 
work and events in the Supernova Cosmology Project (SCP), during the period 1988 to 1998, which 
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PREAMBLE 

In the past I have written papers that described new physics results. In this report, 
however, on this tenth anniversary of the discovery of dark energy, I present the history 
of our discovery as I saw it. I will recount parts of this many-year project that I remember 
particularly vividly as important in the work I was doing. Of course, the other members 
of the team would recall different mixes of events and efforts in which they were 
particularly engaged, and I hope they will find the opportunity to give their recollections. 

When the bubble chamber work in Luis Alvarez's group at LBNL came to an end in 
the early 1970s, there started an interest in astrophysics. Rich Muller became interested 
in the Cosmic Microwave Background (CMB). He was later joined by George Smoot. 
After convincing Luis Alvarez of the feasibility of a CMB measurement and obtaining 
his support, Rich Muller and George Smoot and their collaborators studied the CMB 
with detectors placed on U2 airplane flights. This led to their discovery of the CMB 
dipole asymmetry. Smoot then went on to work with the COBE satellite, and was the 
leader in the discovery of the CMB anisotropy. 

Meanwhile, Rich Muller, joined by Carl Pennypacker (a staff member at the Space 
Science Laboratory), set up an automated search for nearby supernovae (SNe). This 
work was primarily under Pennypacker' s responsibility and was helped by Richard 
Treffers of the Berkeley Astronomy department and others. They were later joined by 
Muller's student Saul Perlmutter who carried out his thesis, a search for "Nemesis", the 
suspected companion star to the sun, as well as developing the analyses and software 
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for the SNe search. During the period 1980-1988, they demonstrated that the method, 
originally suggested by Stirling Colgate et al. [1], worked, and they were able to discover 
about 25 nearby SNe. This work became the prototype of later automated SNe searches. 

In 1988, a National Science Foundation (NSF) Center for Particle Astrophysics 
(CfPA) was being formed on the UC Berkeley campus. Pennypacker and Perlmutter, 
working in Muller's group, had developed a new experiment to discover the "Fate of 
the Universe" through a study of distant Type la SNe, claimed to be "standard candles." 
This experiment was included as one of the elements of the new CfPA [2]. 

THE SCP GROUP AND THE DISCOVERY OF DISTANT SNE 

In 1989, as I was thinking about my next experiment, 2 I was invited by Carl Pennypacker 
to join in the search for the discovery of the "Fate of the Universe" in what was at 
that time called the "Deep Supernova Search." The subject appealed to me, as did the 
proposed technique, since it involved evaluating images, something I have been doing 
throughout my career. 3 

In 1990-91 the entire group at Berkeley (now known as the Supernova Cosmology 
Project (SCP)), consisted of Carl Pennypacker, Saul Perlmutter, Heidi Marvin, myself, 
and Rich Muller (shown in Fig. 1, 1. to r.) . 

Soon after my joining the group, there were major changes in its makeup. In 1991 
Rich Muller decided to spend his efforts in research on the ice ages and global weather 
patterns, and Carl Pennypacker founded the "Hands on Universe" Program for high 
school students and began to devote most of his time to educational activities. This left 
Saul and me to carry on with graduate student Heidi Marvin Newberg and later Alex 
Kim. Shortly thereafter, Bob Cahn, at the time the physics division director at LBNL, 
discussed with me the question of a group leader. With my strong support, he decided to 
appoint Saul to that position. 

Perhaps because Rich, Saul, Carl, and I were not established as astronomers, it was at 
first difficult for us to get time on the large premier telescopes. Before I joined the group, 
on Carl's initiative , Carl, Rich and Saul, together with two colleagues in Australia, Brian 
Boyle (starting out in the UK) and Warrick Couch, did manage to obtain scheduled time 
on the 3.9m Anglo Australian Telescope and built a focus reducing lens and CCD camera 
for it. During the three years we observed at this telescope, while the system worked well 
(and became one of the most-used instruments at the telescope for years after), there was 
no identified SN candidate. Unfortunately, 80% of our scheduled nights were lost due to 
bad weather. The lens and camera construction and installation and the techniques we 
developed, as well as an unidentified candidate found later at the 2.5m Isaac Newton 



2 In 1989 my co-group leader George Trilling and I had come to an end of a nearly twenty-year collab- 
oration with Burton Richter and Martin Perl and coworkers of SLAC and LBL where we had discovered 
the psion family, charmed mesons, and the tau lepton among many others. 

3 Evaluating images is something I have been doing throughout my physics career, beginning with 
photographic emulsions in the 1950s, to bubble chambers in the 1960s, to computer-reconstructed particle 
events in the 1970s, and 1980s. Thus, the supernova experiment with its computer-reconstructed optical 
images, seemed to be a nice fit with my experience and inclinations. 
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FIGURE 1. The SCP group at Berkeley in 1990-91. Carl Pennypacker, Saul Perlmutter, Heidi Marvin, 
Gerson Goldhaber, and Rich Muller 1. to r. 



Telescope, are discussed in great detail in Heidi Marvin Newberg's thesis [3]. 

While for the first three years of our efforts we did not find any identified supernovae, 4 
this time was not wasted, as it allowed us to develop the techniques which led to the 
eventual success of the project. 

In particular, Saul Perlmutter developed the technique for finding "SNe on demand" or 
"batch mode". This involved taking a CCD "reference" image just after a new moon, and 
coming back to take CCD "discovery" images before the next new moon (See Fig. 2). 
The data were sent back to Berkeley the same night they were taken at the telescope. 
On the next morning we ran a program that subtracted the reference image from the 
discovery image and flagged possible candidates. In a few hours of hand scanning we 
were able to select the promising SN candidates. Thus as soon as SN candidates were 
discovered one could send observer(s) to the Keck 10m telescope to measure spectra 
and, after that, start to take data points for the light curve at ground based telescopes. 
(See Fig. 2 for list of follow-up telescopes.) For the first time, new supernovae could 
be guaranteed to be discovered on a certain date, and while they were still brightening- 
making it possible to propose for scheduled nights on the follow-up telescopes. This 



4 However, our general approach received important validation with the discovery by Norgaard-Nielsen 
et al. [4] of a single high z supernova (red shift of 0.31) after a two year effort. This SN was however 
discovered well after maximum light. 
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technique, or some variation of it, was adopted by all subsequent SN searches. 

It is instructive to compare Saul's new technique to what was being invented as the 
state of the art in the early 1990s at much lower redshifts, Hamuy et al.'s Calan/Tololo 
Supernova Search (CTSS). This search observed fields (with photographic plates or 
film) twice a month, and then organized follow-up observation campaigns - but unlike 
Saul's technique is was not aimed at producing guaranteed batches of (on-the-rise) SN 
discoveries on a given date. Thus, towards the end of the CTSS search, Hamuy et al. [5] 
wrote: "Unfortunately, the appearance of a SN is not predictable. As a consequence of 
this we cannot schedule the followup observations a priori, and we generally have to rely 
on someone else's telescope time. This makes the execution of this project somewhat 
difficult." 

As pointed out in a review article by Perlmutter and Schmidt, Saul's technique ad- 
dressed a different problem: "The SCP targeted a much higher redshift range, z > 0.3, 
in order to measure the (presumed) deceleration of the Universe, so it faced a different 
challenge than the CTSS search. The high redshift SNe required discovery, spectroscopic 
confirmation, and photometric follow up on much larger telescopes. This precious tele- 
scope time could neither be borrowed from other visiting observers and staff nor applied 
for in sufficient quantities spread throughout the year to cover all SNe discovered in a 
given search field, and with observations early enough to establish their peak bright- 
ness. Moreover, since the observing time to confirm high redshift SNe was significant 
on the largest telescopes, there was a clear 'chicken and egg' problem: telescope time 
assignment committees would not award follow-up time for a SN discovery that might, 
or might not, happen on a given run (and might, or might not, be well past maximum) 
and, without the follow-up time, it was impossible to demonstrate that high redshift SNe 
were being discovered by the SCP" (Perlmutter & Schmidt [6]). Saul's technique solved 
this problem. 

After we found the first few SNe, Saul was able to convince the various telescope 
scheduling committees to give us time at the matching epochs to carry out this approach. 
In particular this not only allowed us to get more time to do searching, but also to get 
multi-band photometric and spectroscopic follow-up. By 1996 we were even allotted 
time at the Hubble Space Telescope, where, beginning in 1997, we measured lightcurve 
data points for some of the most reliably identified SNe. 

To eliminate cosmic rays and hot pixels, we took two images, usually five minutes 
apart, both for reference images and discovery images. Early in 1990, while scanning 
for SNe candidates (but not finding any!) I did discover a group of about 20 asteroids. 
The asteroids are distinguished by the fact that they are seen to move between the 
two discovery images. However, some of them do not move significantly, and can be 
confused with SNe. We therefore introduced a change in our procedures, interspersing 
a different field between the two discovery images, so as to allow an additional interval 
of 5 to 10 minutes between them. This allowed easy identification of the asteroids. 
Asteroids were later used by Carl Pennypacker in his Hands on Universe program to 
allow students to discover and measure real astronomical objects. 

In 1991 we were joined by visiting scholars Ariel Goobar and Silvia Gabi from 
Sweden and CERN and later on by Reynald Pain from France and Isobel Hook, an 
astronomer with expertise in spectroscopy, from England. 
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Supernova Discovery and Measurement Sequence. 




FIGURE 2. The strategy developed by Saul Perlmutter for finding SNe on demand from repeated CCD 
images. 



In 1993 we collaborated with Alexei Filippenko, of the UC Berkeley Astronomy 
Department, for measurements of spectra at the Keck Telescope in Hawaii. 

In 1994 Don Groom, from the Particle Data Group at LBL, joined our group, as well 
as Susana Deustua. Don Groom adapted the CERN program MINUIT to fit our SNe to 
the light curve and to determine the stretch parameter. 

Graduate students Alex Kim, Matthew Kim and programmer Ivan Small were also 
essential in this early work. 

By 1995, now that we had became successful in finding distant supernovae, we were 
able to hire two postdocs, Rob Knop and Peter Nugent, and later, Greg Aldering. These 
new people, as well as Susana, all had a background in astronomy. 

By 1998 our group had grown to an international collaboration with 32 members 5 



Gregory Aldering, Brian Boyle, Patricia Castro, Warrick Couch, Susana Deustua, Richard Ellis, Se- 
bastien Fabbro, Alexei Filippenko, Andrew Fruchter, Gerson Goldhaber, Ariel Goobar, Donald Groom, 
Isobel Hook, Mike Irwin, Alex Kim, Matthew Kim, Robert Knop, Julia Lee, Chris Lidman, Thomas 
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that signed the Dark Energy discovery paper [7]. In 2007 Saul and the other 31 were co- 
recipients of the Gruber prize in cosmology for the discovery of Dark Energy, together 
with Brian Schmidt and the Hi-Z Team [8]. 

MILESTONES LEADING TO THE DISCOVERY OF DARK 

ENERGY 

Here is a list of some of the important comments and papers which allowed us to make 
our discovery. 

• Suggestions to use Type la SNe as standard candles: From 1930s on Zwicky, Baade, 
Sandage, Kowal, Tammann etc. 

• 1984 Pskovskii [9], 1993 Phillips [10], 1995,1997 Perlmutter et al. [11, 12, 13], 
1995,2001 Goldhaber et al. [11, 12, 14] : light curve shape related to brightness 
and introduction of "stretch". (An alternative approach to this light curve shape 
correction was in 1995, 1996 Riess, Press, and Kirshner [15, 16]. This approach, 
which employed a prior assumption on dust extinction, was used by the Hi-Z team.) 

• 1988 Leibundgut [17]: Type la Lightcurve Template. 

• 1989 Norgaard-Nielsen et al. [4]: First distant SN, z = 0.31 however discovered 
well past maximum light (maximum brightness). 

• 1990 Leibundgut [18], 1993 Hamuy et al. [19]: Single filter K-corrections. To 
convert a filter region of the redshifted spectrum to the same filter at zero redshift. 

• 1992 After three years of searching we discovered and followed our first distant SN 
(z = 0.458) before maximum light [20]. 

• 1993-6 Hamuy, Phillips, Suntzeff, Maza et al. [21, 22]: CALAN/TOLOLO z < 0.1 
reference sample of 29 Type la SNe. Of these we used 18 they discovered before 
10 days past maximum light. 

• 1994-1995 Conference presentations by Perlmutter showing the effectiveness of 
the "SNe on demand" technique. 

• 1995 Kim, Goobar, and Perlmutter [23], 2002 Nugent, Kim, and Perlmutter [24] : 
Cross filter K-corrections. By roughly matching the high-redshift observations' fil- 
ter to the redshifted low-redshift observations' filter this approach solved the prob- 
lem of the large statistical and systematic uncertainites in single-filter K corrections 
due to the large extrapolation errors when the exact form of the spectrum is not 
known and could not be practically measured over the full lightcurve. 

• 1995 Goobar and Perlmutter [25]: Separation of Q. m and Q A shown to be feasible, 
fitting apparent magnitudes of SNe la at a variety of redshifts extending to z > 1 . 
This paper gives the luminosity distance formula for the Hubble plot calculations 



Matheson, Richard McMahon, Heidi Newberg, Peter Nugent, Nelson Nunes, Reynald Pain, Nino Pana- 
gia, Carl Pennypacker, Saul Perlmutter, Robert Quimby, Pilar Ruiz-Lapuente, Brad Schaefer, Nicholas 
Walton. (Alexei Filippenko moved to the Hi-Z Team when it was formed in 1995.) 
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and shows that at each redshift the best-fit confidence region will make a diagonal 
on the CL m vs. £l A plane (with a slope that varies with redshift). 

Here "stretch", introduced by Saul, is the deviation of the measured timescale of a 
lightcurve from a standard average lightcurve. Typical values are s = 0.8 to 1.2. The 
larger s the brighter the SN. The effective B-magnitudes are corrected to a standard 
lightcurve. 

SEPT. 24, 1997: A PEAK IN THE £l m HISTOGRAM 

What has become the conventional way to determine Q. m and Q.^ is the Goobar- 
Perlmutter [25] best fit confidence level distribution on the Q. m , plane. This approach, 
favored by astronomers, requires accurate determinations of the errors and correlated er- 
rors. When the SNe light curves are first fitted to a template this information is not yet 
available. 

The alternative approach, similar to the technique used to find resonances in particle 
physics, is to look for peaks in the distributions (on histograms) of the variable being 
studied, namely mass in particle physics and Q. m in our case here. In this approach the 
measurement errors are reflected in the width of the resulting peak. 

In September 1997 we had completed a first pass of the light curve data point analysis 
on 38 SNe, but the detailed error analyses were not yet available. It took several months 
from the time a SN was discovered until the light curve points could be evaluated. We 
also took some final reference points a year after discovery. Rob Knop, who was working 
on this, was able typically to complete the measurements on one to two SNe per month. 

On the basis of these 38 Type la SNe we obtained an indication that the universe 
is accelerating, rather than decelerating as we had originally expected. This was the 
culmination of nine years of work learning how to find high redshift SNe, how to 
measure them, how to fit them, how to K-correct them, how to stretch correct them, 
how to account for dust reddening, and then finally fit them to a lightcurve template. 
This was all achieved by a collaborative effort of the entire group. I studied the SNe 
after the data points on the lightcurve had been measured, fitted them with the group's 
lightcurve fitter which gave us the stretch value, and made a table with some 20 attributes 
for each SN. 

The question was next how to obtain the Q. m distribution from this plot. What I tried 
was to plot Hubble curves, for a flat universe, with Q. m in fine intervals at 0.0, 0.2, 0.4, 
0.6, 0.8, 1.0. See Fig. 3. A distribution in Q. m is then obtained by counting (by hand!) 
all the SNe that fell into each Ml m interval giving the histogram shown in Fig. 4. A few 
days later Saul wrote a program, SNLOOK, to have the computer count the number of 
SNe in each AQ. m interval. He also made plots both for a flat universe, and for a universe 
with Q. A = 0. 

Here it should be noted that a flat universe was not completely established at that point 
in time (1997). (The flat universe case was based in part on inflation theory [26, 27, 28] 
and in part on very early rumors from as yet unpublished CMB fluctuation data [29] . The 
detailed measurements of the fluctuations in the CMB, demonstrating a flat universe, 
came later [30, 31, 32].) However, an analysis of the supernova data for the cases of a 
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FIGURE 3. Our preliminary data of 38 SNe, as of Sep., 24, 1997. The effective B-magnitude vs. log(cz). 
Here the effective B-magnitude is the K-corrected and stretch-corrected magnitude and z is the redshift. 
The curves for a series of £l m values, in a flat universe, are shown. The sum of the data points in each Ai2,„ 
interval, as counted by hand, is also shown at the right hand upper edge of the figure. 



flat universe and a Qa = universe is a very informative way to characterize the data 
and its implications. As the 1995 Goobar and Perlmutter paper[25] had shown , the 
measurements of high-redshift supernovae in the redshift range we were studying result 
in a long diagonally-oriented confidence region in the Q. m — Q A plot, like the one shown 
in Figure 9. We can also analyze this data by projecting onto an Q. m histogram for either 
of the simple cosmology cases, the flat universe line or the Qa = line in Figure 9. Given 
the ~45-degree slope of the diagonal confidence region for the SN redshift ranges we 
were studying, if we find a very low value of Q. m in the flat-universe-case histogram, 
we will be guaranteed to find a negative value for Q. m in the ^a = 0-case histogram. 
Since there is only positive mass in our physics, this implies that there is a cosmological 
constant or some other form of Dark Energy - and the result holds whether or not the 
universe is actually flat! This is very different from the low Q. m measurements from 
galaxy clusters, etc., that were being made in the late 1990s, where the measurement by 
itself did not tell you anything about 

In Fig. 4 the peak we actually saw at Q. m = 0.2 to 0.3 for a flat universe indicates 
that Qa is 0.8 to 0.7. Thus we obtained a significant positive value for Qa or possibly 
a cosmological constant. As we were looking for deceleration we always calculated the 
deceleration parameter qo = ( l/2)^ m — ^a- Here go = —0.7 to —0.55. A negative value 
for qo implies acceleration. We calculated go for each of our samples and it always said: 
ACCELERATION. 
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FIGURE 4. The £l m distribution for the data points in Fig. 3, as presented at our SCP group meeting on 
Sep. 24, 1997, is shown. The values for the first 7 SNe, which gave a considerably larger value for £l m , 
are shown cross-hatched. Here £1a is given by 1 — Cl m . 



The low redshift distribution of SN absolute magnitudes available at that time already 
showed a strong peak that indicated a large fraction of supernovae were not suffering 
extinction. The strong peak in the Q. m distribution was likely also to be due to unex- 
tincted SNe, but we understood that, whatever this extinction distribution, as long as 
it was the same at low and high redshift the cosmology results would be correct. We 
studied the color distributions at both redshift ranges, and found that the distributions 
were sufficiently consistent. This histogram of Figure 4 (and the others shown below ) 
therefore accounts for dust in this statistical approach, rather than correcting extinction 
SN by SN. (Our other analyses at that time examined this point with further studies of 
the measured colors of the distant SNe, performing analyses correcting individual SN 
brightnesses and analyses comparing groups of SNe. We concluded that our cosmology 
result was robust - and that the correct way to analyze color and dust for all the data sets 
at that time was not to use a Baysean prior on the dust distribution, as others had done.) 
It is also interesting to note that the stretch correction, which becomes very important in 
later precision measurements, does not change these preliminary values significantly. It 
only tends to broaden the observed peak. At any rate we did apply the stretch correction 
to these data. 
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This startling result was naturally treated with some skepticism within the group, 
when these plots were discussed in the weekly SCP meetings at the end of September 
1997. Although we had already submitted a Nature paper[33] indicating evidence for 
a cosmological constant, we were still surprised that this result was settling in on the 
lower end of the Q. m range from our first seven SNe [11, 12, 13] which had centered on 
a high value for Q. m . 6 These seven SNe are shown as cross hatched in Fig. 4. and indeed 
occur at the upper end of the distribution. 

Again the fact that the deceleration coefficient was negative, indicating acceleration, 
was hard to swallow at first. I had been "bump hunting" for the past 30 years, and found 
the observed peak completely convincing. I believe that Saul convinced himself after 
he wrote the program, SNLOOK and by independently studying our data. He later in 
October and December presented such a histogram in colloquia giving Q. m and £1^. 
Some other colleagues, Rob Knop and Carl Pennypacker, appeared convinced; while 
Greg Aldering was not fully convinced, he stated however that this histogram "helped 
to galvanize the effort within our group". Indeed, as the result of a major effort, the SCP 
group was able to calculate the best fit confidence levels, for all 40 SNe, in time for the 
Jan. 8, 1998 AAS meeting. 

For confirmation, I asked my colleagues to check these results, in case there could 
have been a mistake. However, all the 20 entries in the so called "Gerson-table" for each 
of the 38 SNe were confirmed as correct. Later Greg Aldering added more columns to 
this table, giving the correlated errors, and Richard Ellis added a column with the best 
estimate of the nature of each host galaxy. 

Although as mentioned above, Goobar and Perlmutter (1995) [25] had shown that any 
SN measurements indicating low Q. m for a flat universe will also indicate negative £l m 
for a £2a = universe, it was striking to see this directly in the data. Thus, in Fig. 5. the 
Hubble curves for both a flat universe and an Qa = universe are shown superimposed. 
For an Q. A = universe we obtain a corresponding peak on the Q. m histogram. However 
now it occurs for negative Q. m values. In Fig. 6. is shown the Cl m histogram, this time 
calculated with Saul's SNLOOK program, for this case. As can be noted, 25 out of the 
38 SNe give a negative and hence un-physical value for Q. m . This clearly demonstrates 
that Qa cannot be zero but must be greater than zero. Figures 5 and 6 are from my 
notebook from October 1997 and were shown informally to group members. This result 
was also presented in my seminar at Santa Barbara in December, 1997. 

One can also turn this argument around and calculate the location of the peak in an 
Q. m histogram for a series of different increasing Qa values and find for which Qa value 
the peak begins to lie at positive (hence physical) values. This approaches, and goes 
beyond, £l m , Q.^ values consistent with the above values for a flat universe. This method 
is equivalent to traveling along the central line for the confidence level "ellipses", starting 
at negative Q. m values, shown in Fig. 9. 



In retrospect, it is instructive to revisit the question of what happened with those first seven supernovae. 
As shown shaded in Fig. 4 they are all at high Q.,„ values. Two turned out to be outliers with one probably 
not a Type la. The other 5 moved slightly towards lower mass densities with remeasurement, but for the 
most part they simply happened to lie on the tail of the distribution. The lesson is: beware of statistics of 
small numbers! 
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FIGURE 5. The Hubble curves for both a flat universe (solid curves) and an Q.^ = universe (dotted). 




FIGURE 6. The Q. m distribution for a H A = universe. Note that 25 of the 38 SNe have a negative Q. m 
and hence an un-physical value. This shows clearly and independently of the flat universe condition that 
il A must be greater than zero. 
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TABLE 1. SCP SNe la Discoveries with z > 0.2 



Epoch 


Number 


Total 


Telescope 


1989-91 








Anglo Australian 3.9m 


1992 


1 


1 


La Palma 2.5m 


1993 


6 


7 


Kit Peak 4m 


1995 


9 


16 


Cerro Tololo 4m 


1996 


7 


23 


Cerro Tololo 4m 


1997 


19 


42 


Cerro Tololo 4m 



We thus have two ways to visualize with histograms the evidence for a positive value 
for The first histogram essentially looks at where the diagonal confidence region (in, 
e.g., Fig. 9) crosses the flat universe line, and the second looks at where it crosses the 
£1\ = axis. Both indicated that Q A must be positive (since the £l A = histogram gives 
unphysical results). And, of course, when we calculated qo the supernova data clearly 
indicated acceleration. 

In Table 1 are shown the SN discoveries over the years 1989 to 1997. Saul's SNe on 
demand method really came into its own when we obtained time on the Cerro/Tololo 4m 
telescope combined with the "big throughput" camera (BTC) of Bernstein and Tyson. 
Table 1 shows that over a nine-year period we discovered 42 SNe with z > 0.2. 

Oct. 1997. Two papers by the SCP and Hi-Z Team as well as a joint 
paper submitted with hints, from very small samples of SNe, for a 

non-zero Q.a 

The very first analysis of multiple high-redshift SNe indicated a high Q. m value, 
albeit with large statistical error bars: using the 5 of these first 7 SNe with lightcurve 
timescales that were calibrated at low redshift yielded Q. m = 0.94+Q2g for a flat universe 
[11, 12, 13]. The following year, we added one very-well-observed type la supernovae 
that we discovered and followed with the Hubble Space Telescope (HST), SN 1997ap, 
that was at that time again the most distant that had ever been observed, z = 0.83, with 
the most leverage for this measurement. The addition of this single SN allowed us to 
recalculate £l m for 5+1 supernovae, giving Q. m = 0.6 ±0.2 for a flat universe - for the 
first time in a publication a direct indication of a positive cosmological constant, albeit 
only at the two sigma level. Fig. 7 shows the Hubble diagram from this paper, where the 
new well-measured highest redshift SNe indicates a positive cosmological constant with 
even a stronger significance. We submitted this paper to Nature on Oct. 6, 1997 [33] and 
it appeared in the January 1, 1998 issue. 

Garnavich et al. [34] submitted a paper by the Hi-Z Team on Oct. 13, 1997, based 
on three confirmed high redshift SNe measured with the HST, and one measured from 
the ground. They concluded that matter alone is insufficient to produce a flat universe. 
Specifically for Q. m + Qa = 1 > is less than 1 with more than 95% confidence, and 
their best estimate of Cl m is —0. 1 ± 0.5 if ^a = 0. 

A joint paper by some members of the two groups Riess et al. [35] submitted Nov. 
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21, 1997 and accepted Apr. 17, 1998, showed that the "snapshot" method works. That 
involves taking a single photometry point and a single spectrum on the same night and 
deducing the magnitude at maximum light from this limited information. On the basis 
of four SNe handled in this fashion, as well as the SNe from the above two papers, they 
deduced that for a flat universe Q m = 0.19+q^ and Q. m = — OJI+q^ for an £l A = 
universe. 

While these SCP and Hi-Z papers, published in 1998, gave hints of a non-zero Q.^ 
they relied on a very small number of SNe. The hints from the snapshot paper similarly 
relied on four SNe using this less precise method. 

All this work was going on more or less contemporaneously with my study of the 
38 SNe. My colleague Peter Nugent, who did some of the fits of these confidence 
level distributions, felt that these 3 "hints" were a more convincing and possibly earlier 
evidence for a non-zero than the peak in the Q. m histogram. My feeling is that we 
needed the larger statistics since we had amply demonstrated that with a very small 
number of SNe (7 SNe) one can be way off! 



OCT. 23 TO DEC. 14, 1997: THREE COLLOQUIA AND A 

SEMINAR 

The first four public presentations of the SCP evidence for Q.^ > and acceleration, 
based on the peak in the Q. m histogram for 40 SNe as well as the unphysical results 
obtained when Q.^ is assumed to be zero. 



13 



From October to November, we re-measured and re-fitted all of the 38 SNe and added 
two more, giving 40 (while three were identified as outliers). We also changed from 
the Burstein and Heiles [36] version of dust extinction in our galaxy to the more recent 
Schlegel, Finkbeiner and Davis [37] version. I then revised the tables to reflect all these 
changes. Both Saul and I plotted the Hubble plots, based on the revised tables, and used 
SNLOOK to obtain the latest version of the histogram. 

The first public presentation of our results was a colloquium by Saul on October 23, 
1997 at the Physics Department, UC San Diego. This was followed by Saul's colloquia 
at the Physics Department UC Berkeley on Dec. 1, 1997 and at the Physics Department 
UC Santa Cruz on Dec. 11, 1997 (See Fig. 8). 



Supernova Cosmology Project 
Corrected for 




Qm (flat universe) 

FIGURE 8. {left) The figure of the £2„, distribution shown at Colloquia at UC San Diego in October 
1997 and at UC Berkeley and Santa Cruz in early December 1997 by Saul Perlmutter. The x-axis gives 
Q. m for a flat universe. All the SNe given in Fig. 4. have been re-fitted. The Gaussian fit to this preliminary 
data gave £2„, = 0.29 ±0.06 with a a = 0.3. This analysis accounts for dust with the more correct statistical 
approach discussed in the text above concerning Fig. 4, rather than correcting extinction SN by SN. (right) 
A reproduction of the actual figure of the £2„, distribution shown at the ITP Santa Barbara on Dec. 14, 1997 
by Gerson Goldhaber. The lower scale gives £l m for a flat universe. The upper scale gives £2„, for an Q.^ = 
universe. Note that this assumption gives negative (and thus unphysical) values for £l m . This demonstrates 
that i^A must be greater than zero by a method which does not require the flat universe condition. 

Next, I gave a seminar at the Institute for Theoretical Physics (ITP) at UC Santa 
Barbara on Dec. 14, 1997 (See Fig. 8). In all four of these talks, as stated above, the 
conventional description of the confidence level distribution on the Q. m , plane was 
not yet available for the 40 SNe. So in all four talks the Q, m histograms, re-evaluated 
with more careful re-measurements, were shown. We showed that the result was robust 
against the likely sources of systematic error, but we still considered our result as 
preliminary as we had not yet completely explored every possible effect. 
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It was not easy to convince the astrophysics community of this result, as it was 
contrary to all ingrained beliefs. 7 Acceleration rather than deceleration as expected! 

The question has been posed: who remembers the details of these talks? I know that 
Kim Griest remembers Saul's talk at UC San Diego and was enthusiastic about the result. 
I also know Joel Primack and Michael Riordan remember Saul's talk at UC Santa Cruz. 
As former and current particle physicists they understood the significance of a peak on 
a histogram. Joel Primack was particularly enthusiastic and stressed the importance of 
this discovery after Saul's colloquium. Dave Branch remembers my talk and subsequent 
discussions and understood its significance. David Gross asked me after my talk how I 
could come to such a momentous conclusion on the basis of just 40 SNe. 

As is usual, Saul's Colloquium at UC Berkeley was video taped (DVD available). 
Saul first talked about the "hint" from our 5+1 SNe fit I mentioned above. Then based 
on a gaussian fit to the observed peak on the histogram, Saul quoted Q. m = 0.29 ± 0.06 
for a flat universe and thus £l A = 0.71. In my presentation at the ITP I did not try for an 
accurate value of Cl m , but rather showed Fig. 8 and indicated that Q. m = 0.3 agrees with 
our data for a flat universe. In addition I demonstrated what happens when we assume 
= 0. As seen from the upper scale in Fig. 8 this assumption gives negative and hence 
unphysical values for Q. m . As stated earlier, given the slope of the confidence region for 
any supernova data, such values for Q. m in the two cases considered demonstrate that Q.^ 
must be greater than zero, and does not depend on assuming a flat universe. 

In his talk Saul went on to point out that while on the one hand 0.7 for Q.^ corre- 
sponded to a large fraction of the universe but on the other hand it was a very small 
value, by a factor 10 122 , compared to the value expected from virtual particle vacuum 
fluctuations. 

It is remakable that in Saul's talk on Dec. 1, 1997 we had already presented, in public, 
evidence for the whole story of what was later called Dark Energy by Michael Turner. At 
DM08 I presented a four-minute video excerpt of Saul's talk, showing him making these 
points based on the histogram. (This video is given in the DM08 conference proceedings 
[41].) 

Not much has changed qualitatively in the 10 years since then. With many more SNe 
studies [42] as well as concordance with CMB and evidence from cluster studies [43] 
and later Baryon Acoustic Oscillation (BAO) [44] studies, we now know Q. m to more 
significant figures, but qualitatively there is no change from the value we quoted in 
1997. We have learned that the equation of state parameter w is very close to -1 and have 
some limits on the variation of w with z. But after all that, just as 10 years ago, we still 
do not understand the nature of Dark Energy [45] . 



Examples of exceptions to the il A = assumption are: Steve Weinberg [38] made an argument, on the 
basis of the anthropic principle, that i2 A could be comparable to the currently observed value. Krauss and 
Turner [39] argued primarily on the basis of the age of the universe and other then available cosmological 
measurements that i2 A should lie between 0.6 and 0.7. Ostriker and Steinhardt [40] derived a concordance 
model based on all the then available cosmological data and concluded i2 A = 0.65 ± 0. 1 . 
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JAN. 8, 1998: FIRST PUBLIC PRESENTATION OF THE SCP 
RESULTS AT AN AMERICAN ASTRONOMICAL SOCIETY 

MEETING 

The fifth public presentation of the SCP evidence for Q.^ > and acceleration, based on 
the fit to the confidence level distribution on the Q. m , plane for 40 SNe. 

From October to late in December 1997, the SCP group at Berkeley re-measured and 
re-checked and completed all error and correlated error calculations on what had by 
then become 40 SNe. This effort was led by Saul, Greg Aldering, Peter Nugent, and 
Rob Knop. Each group member as well as the group members in Baltimore, Stockholm, 
Paris, England, Australia, Spain and Chile convinced himself/herself by direct study 
of the data that indeed Q. A > and that contrary to preconceived notions we were 
dealing with acceleration of the universe, namely a negative deceleration coefficient. 
The collaboration was prepared to make a major presentation. 

The following month, on Jan. 8, 1998 [46, 47], Saul Perlmutter presented our data 
at the A AS meeting in Washington DC. By this time we had the much more elaborate 
analysis, based on detailed error calculations, by the entire SCP group. We had 40 fully 
analyzed SNe and the confidence level calculation was already available on the Q. m , 
Q.A plane (Fig. 9). Here it should be noted that the confidence level distribution does 
not require the flat universe condition to give Q.a greater than zero. This was the first 
presentation at an official conference of evidence for Q.^ > 0, and hence from finding 
qo < 0, the acceleration of the expansion of the universe. To take possible systematic 
uncertainties into account, we allowed for a very generous systematic error and they 
were shown for the case where all systematics conspired in the same direction. This is 
shown by the dotted confidence level curves in Fig. 9. This figure was reproduced by 
James Glanz in his report on our presentation in Science [48]. Furthermore, using the 
same Q. m histogram for a flat universe as in Fig. 8, Saul also showed the studies of the 
effects of systematics: Malmquist bias, stretch dependence, and color/dust dependence. 

FEB. 18, 1998: PUBLIC PRESENTATIONS AT THE DARK 
MATTER 1998 (DM98) CONFERENCE. 

The sixth and seventh public presentations of the SCP evidence for Q.^ > and accel- 
eration based on 42 high redshift SNe were given. The first public presentation by the 
Hi-Z Team of evidence for Q.\ > and acceleration based on 14 independent high red- 
shift SNe was also given. Both teams quoted results based on fits to the confidence level 
distribution on the Q m , Q.\ plane. 

The following month, on Feb. 18, 1998, both I and Saul Perlmutter gave talks, in that 
order [49], showing our results at the meeting organized by Dave Cline just 10 years ago 
(Dark Matter 1998) at Marina del Rey. By this time our number of fully analyzed SNe 
had grown to 42. ( See Table 1.) 

I spent part of my talk on some work I had done on time dilation in SNe explosions, 
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No Big Bang 




FIGURE 9. The £2a vs Q m best fit confidence level distribution shown at the Washington meeting of 
the AAS on Jan, 8 1998 by Saul Perlmutter. The dashed curves represent a very generous estimate of what 
possible systematic errors could do. 



proving that the redshift is due to the expansion of the universe, rather then to the "tired 
light" hypothesis [11, 12, 14]. I then showed the best fit confidence level distribution on 
the Q. m , Q\ plane, which for a flat universe implied Q. m = 0.281qq8 (statistical) Ijj 04 
(systematic) and hence, from evaluating qo, acceleration of the universe at the present 
epoch. 

In his talk Saul Perlmutter gave more details on our results and on the error calcula- 
tions, and stressed the studies of the systematic uncertainties. 

Following our talks, Alexei Filippenko presented the results [50] of the Hi-Z Team, 
who claimed they had established the acceleration of the universe on the basis of a 
confidence level analysis of 16 high redshift SNe. These consisted of 10 well measured 
SNe, plus 4 from the "snapshot" method [35] plus 2 which came from our set of 42 [7]. 

Later the refereed publications of the two groups capped these momentous results. 
The Hi-Z Team paper was submitted on Mar. 13, 1998 [8] and the SCP paper was 
submitted on Sept. 8, 1998 [7]. 
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HOW DOES AN £l m HISTOGRAM LOOK TODAY? 



To see how the method, patterned on particle physics and used to indicate acceleration 
on Sept. 24, 1997, works with larger statistics I have analyzed (this time with a computer 
program) a sample of 257 SNe with z > 0.2 [28] in the same fashion as in Figs. 3 and 4. 
Fig. 10 shows the resulting distribution in Q, m . The enormous peak is clearly observed, 
which completely confirms the validity of the peak for the 38 SNe in Fig. 4. 

Omega_M from Mag at Max S & C CORRECTED 

80 1 | i i i i | i i i i | i i i i | i i i i | ' ' ' 
Chi'iq = 0.70 Om_M = 0.29 Sigma = 0.29 Max cnt : 



8.875 




0.0 0.5 1.0 

OmegQ_M Flat Universe 



FIGURE 10. The D.„, distribution for a sample of 257 SNe with z > 0.2 out of a total of 307 which were 
collected and re-analyzed in our group by Marek Kowalski, David Rubin et al. [29] in connection with 
a current (2008) SCP paper giving a Union of all published SNe. The SN magnitudes at maximum light 
were K-corrected as well as stretch and color corrected. Here as above £2a is given by 1 — £l m . 



FEB. 20, 2008: SO WHAT IS THE "FATE OF THE UNIVERSE" ? 

By now, 10 years later, there are many experiments which have confirmed Dark Energy. 
We can however still ask: what is the "Fate of the Universe" ? 

The question is now in the form: is constant or does it vary with z? Or alternatively: 
is w = — 1, namely are we dealing with Einstein's cosmological constant ? Or can one 
measure deviations from — 1 where w is the equation of state parameter? 

We hope the Joint Dark Energy Mission (JDEM) will answer this question before the 
next decade is up at DM2018. 
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SUMMARY OF THE STEPS IN THE SCP DISCOVERY OF 
EVIDENCE FOR DARK ENERGY 



• 1988 first pilot search for high-redshift SNe at Palomar organized by Carl Penny- 
packer and Saul Perlmutter. 

• 1988 proposal to NSF center at UCB for "Fate of Universe" study. 

• 1992 discovery of SN1992bi, z = 0.486. 

• 1995 first 7 SNe yielded large Q. m value. 

• Sept. 24, 1997 Q. m = 0.2 to 0.3 for flat universe from peak in histogram, SCP group 
meeting based on 38 SNe. 

• Oct. 1997 hints for non zero Qa from SCP and Hi-Z Team. 

• Oct. 23 to Dec. 14, 1997 public talks by Saul Perlmutter and Gerson Goldhaber 
showing evidence that for a flat Universe Q. m = 0.3 and hence ^a = 0.7 based on 
the peak in the Q. m histogram, Figs. 8, and by then on 40 SNe. 

• Jan. 8, 1998 SCP presentation by Saul Perlmutter at the AAS meeting in Wash- 
ington DC showing the Q. m , Q./^ plane with best fit confidence level distributions 
yielding flat universe values of Q. m = 0.3 and Qa = 0.7 . 

• Feb. 18, 1998 Both groups talk at Marina del Rey DM98. Gerson and Saul show 
evidence for a non-zero ^a and acceleration based on 42 SNe, followed by Filip- 
penko who showed evidence for a non-zero Q.^ and acceleration based on 10 + 4 + 
2 SNe for the Hi-Z Team. 

ACKNOWLEDGEMENTS 

I wish to thank Kyle Barbary, Alex Gude, David Rubin, Tony Spadafora and Nao Suzuki 
for help with this manuscript. 

This work has been supported in part by the Director, Office of Science, Office of 
High Energy Physics, of the U.S. Department of Energy under Contract No. DE-AC02- 
05CH11231. 

REFERENCES 

1. S. A. Colgate, E. P. Moore, and R. Carlson, PASP 87, 565 (1975). 

2. C. Pennypacker, and R. A. Muller, Measurement of the gravitational mass density (1988), a Proposal 
to the National Science Foundation. 

3. H. J. M. Newberg, Measuring qO using supernovae at Z approximately 0.2, Ph.D. thesis, California 
Univ., Berkeley. Lawrence Berkeley Lab. (1992). 

4. H. U. Norgaard-Nielsen, L. Hansen, H. E. Jorgensen, A. Aragon Salamanca, and R. S. Ellis, Nature 
339,523(1989). 

5. M. Hamuy, et al., AJ 106, 2392 (1993). 

6. S. Perlmutter, and B. P. Schmidt, "Measuring Cosmology with Supernovae," in Supernovae and 
Gamma-Ray Bursters, edited by K. Weiler, 2003, vol. 598 of Lecture Notes in Physics, Berlin 
Springer Verlag, p. 195. 

7. S. Perlmutter, et al., ApJ 517, 565 (1999). 

8. A. G. Riess, et al., AJ 116, 1009 (1998). 



19 



9. Y. P. Pskovskii, Soviet Astronomy 28, 658 (1984). 

10. M. M. Phillips, ApJL 413, L105 (1993). 

11. S. Perlmutter, et al., and G. Goldhaber, et al. Presented at NATO Advanced Study Inst, on Thermonu- 
clear Supernovae, Aiguablava, Spain, 20 - 30 June 1995. 

12. P. Ruiz-Lapuente, R. Canal, and J. Isern, "Thermonuclear supernovae. Proceedings.," in NATO ASIC 
Proc. 486: Thermonuclear Supernovae, 1997. 

13. S. Perlmutter, et al., ApJ 483, 565 (1997). 

14. G. Goldhaber, et al., ApJ 558, 359 (2001). 

15. A. G. Riess, W. H. Press, and R. P. Kirshner, ApJL 438, 17 (1995). 

16. A. G. Riess, W. H. Press, and R. P. Kirshner, ApJ 473, 88 (1996). 

17. B. Leibundgut, Light curves of supernovae type, I., Ph.D. thesis, Univ. Basel. 137 (1988). 

18. B. Leibundgut, A&A 229, 1 (1990). 

19. M. Hamuy, M. M. Phillips, L. A. Wells, and J. Maza, PASP 105, 787 (1993). 

20. S. Perlmutter, et al., ApJL 440, L41 (1995). 

21. M. Hamuy, et al., AJ 109, 1 (1995). 

22. M. Hamuy, et al., AJ 112, 2391 (1996). 

23. A. Kim, A. Goobar, and S. Perlmutter, PASP 108, 190 (1996). 

24. P. Nugent, A. Kim, and S. Perlmutter, PASP 114, 803 (2002). 

25. A. Goobar, and S. Perlmutter, ApJ 450, 14 (1995). 

26. A. H. Guth, Physical Review D 23, 347 (1981). 

27. A. D. Linde, Physics Letters B 108, 389 (1982). 

28. A. Albrecht, and P. J. Steinhardt, Physical Review Letters 48, 1220 (1982). 

29. E. Torbet, et al., ApJL 521, L79 (1999). 

30. S. Hanany, et al., ApJL 545, L5 (2000). 

31. P. de Bernardis, et al., Nature 404, 955 (2000). 

32. C. L. Bennett, et al., ApJS 148, 97 (2003). 

33. S. Perlmutter, et al., Nature 391, 51 (1998). 

34. P. M. Garnavich, et al., ApJL 493, L53 (1998). 

35. A. G. Riess, P. Nugent, A. V. Filippenko, R. P. Kirshner, and S. Perlmutter, ApJ 504, 935 (1998). 

36. D. Burstein, and C. Heiles, AJ 87, 1 165 (1982). 

37. D. J. Schlegel, D. P. Finkbeiner, and M. Davis, ApJ 500, 525 (1998). 

38. S. Weinberg, Reviews of Modern Physics 61, 1 (1989). 

39. L. M. Krauss, and M. S. Turner, General Relativity and Gravitation 27, 1 137 (1995). 

40. J. P. Ostriker, and P. J. Steinhardt, Nature 377, 600 (1995). 

41. S. Perlmutter, Colloquium (1997), presented at the Universtiy of California, Berkeley. A 5 minute 
excerpt is avalible at http://www.physics.ucla.edu/hep/hep/dm08/talks.html. 

42. M. Kowalski, et al., ApJ 686, 749 (2008). 

43. N. A. Bahcall, R. Cen, R. Dave, J. P. Ostriker, and Q. Yu, ApJ 541, 1 (2000). 

44. D. J. Eisenstein, et al., ApJ 633, 560 (2005). 

45. D. Rubin, et al., ArXiv e-prints (2008). 

46. S. Perlmutter, et al., ArXiv Astrophysics e-prints (1998), poster presented at the American Astronom- 
ical Society meeting, Washington, D.C., January 8, 1998. 

47. S. Perlmutter, et al., "Cosmology From Type IA Supernovae: Measurements, Calibration Techniques, 
and Implications," in Bulletin of the American Astronomical Society, 1997, vol. 29, p. 1351. 

48. J. Glanz, Science 279, 651 (1998). 

49. G. Goldhaber, and S. Perlmutter, Physics Reports 307, 325 (1998). 

50. A. V. Filippenko, and A. G. Riess, Physics Reports 307, 31 (1998). 



20 



